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Abstract We present an example of a boundary layer tailward of the dawn terminator which is entirely 
populated by rolled-up flow vortices. Observations were made by Wind on 24 October 2001 as the 
spacecraft moved across the region at X — 13 R E . Interplanetary conditions were steady with a near-radial 
interplanetary magnetic field (IMF). Approximately 1 5 vortices were observed over the 1 .5 h duration of 
Wind's crossing, each lasting ~5 min. The rolling up is inferred from the presence of a hot tenuous plasma 
being accelerated to speeds higher than in the adjoining magnetosheath, a circumstance which has been 
shown to be a reliable signature of this in single-spacecraft observations. A blob of cold dense plasma 
was entrained in each vortex, at whose leading edge abrupt polarity changes of field and velocity 
components at current sheets were regularly observed. In the frame of the average boundary layer velocity, 
the dense blobs were moving predominantly sunward and their scale size along X was ~7.4 R E . Inquiring 
into the generation mechanism of the vortices, we analyze the stability of the boundary layer to sheared 
flows using compressible magnetohydrodynamic Kelvin-Helmholtz theory with continuous profiles for the 
physical quantities. We input parameters from (i) the exact theory of magnetosheath flow under aligned 
solar wind field and flow vectors near the terminator and (ii) the Wind data. It is shown that the 
configuration is indeed Kelvin-Helmholtz (KH) unstable. This is the first reported example of KH-unstable 
waves at the magnetopause under a radial IMF. 


1. Introduction 

There is a long history of observations of waves at the boundary between the magnetosphere and the mag- 
netosheath [e.g., Lepping and Burlaga, 1 979; Sckopke et al., 1 981 ; Chen and Kivelson, 1 993; Farrugia et al, 
2001 , and references therein]. In view of the velocity shear that exists between these two plasma regimes, 
the Kelvin-Helmholtz (KH) instability has often been invoked to explain these waves. 

From theoretical studies of the KH instability, two main points to keep in mind are the following: (i) the mag- 
netic tension force (analogous to the surface tension force in hydrodynamics) and (ii) the compressibility of 
the plasma. Both are stabilizing factors. Thus, KH instability depends on the magnetic field configurations, 
in particular their orientations with respect to the flow, and the speed of the plasma, which increases with 
distance down the flanks. 

KH waves are thought to be one way of transferring solar wind momentum and energy to the magneto- 
sphere. The KH instability forms part of the so-called "viscous-type" solar wind-magnetosphere interactions, 
to distinguish them from reconnection between the magnetosheath and magnetosphere fields. The con- 
tribution of viscous-type interactions to the cross-polar cap potential is often estimated as ~30 kV [Cowley, 

1 982]. The question of magnetosheath mass entry goes beyond considerations of ideal MHD stability since 
other processes are required to break the associated frozen-in condition. However, the large vortices gen- 
erated by the KH instability may set up conditions favorable to small-scale tearing of magnetic field lines 
inside the structures and, as a consequence, to mass diffusion [see, e.g., Otto and Fairfield, 2000; Smets et al., 
2002; Otto and Nykyri, 2003]. 

Many of the data examples of KH instability in a magnetospheric context have used the capability of 
multiple spacecraft observations, such as Cluster, to confirm the presence of waves and their features, in 
particular if they have reached a nonlinear phase and started to roll over. However, we do not often have 
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this luxury and there are many 
tabulated crossings of a wavy mag- 
netopause boundary made by single 
spacecraft. Can we somehow infer the 
presence of rolled-up vortices from 
single-spacecraft observations? 

A key advance in this direction was 
made by Takagi et al. [2006], Their 
MHD simulations showed that in situ 
observations of a low-density magne- 
tospheric plasma moving tailward at 
speeds higher than that of the adjacent 
magnetosheath is a very good indi- 
cator of rolled-up vortices. This opens 
new possibilities. First to apply this cri- 
terion were Hasegawa et al. [2006], who 
confirmed results on rolled-up KH vor- 
tices obtained earlier by Hasegawa et 
al. [2004] with a multispacecraft anal- 
ysis. While the simulations were done 
fora northward interplanetary mag- 
netic field (IMF) and specific parameters 
characterizing the ambient regions, 
Nakamura et al. [2004] had already 
given the physical origin of the sig- 
nature of a rolled-up vortex. For force 
balance to hold across the vortex (same 
centrifugal forces at a given radial distance from the center), the hot tenuous plasma must revolve at a 
higher speed than the cold dense plasma. 

Figure 1 presents a schematic to help visualize this point. Figure 1 (top) illustrates the perturbed magne- 
topause (MP) at the equatorial dawn flank that begins to roll over into a vortex by the KFH instability. The 
magnetosheath flow is tailward [V x < 0), while the magnetosphere is stagnant. Accordingly, across the 
boundary layer there is a velocity gradient. The drawing is shown in the frame of the vortex, so that the cold 
dense magnetosheath tongue (blue) protruding to the left is slowing down relative to the average flow, 
while the related hot tenuous magnetosphere protuberance (red) points to the right and accelerates toward 
the tail. The cold dense plasma intermingles with the hot tenuous plasma. The thick arrowed lines give an 
indication of the plasma motion. The thin blue lines are the conjectured deformation of the magnetosheath 
magnetic field projected into the XT plane. Flence, (1) we expect an alternation of high- and low-density 
cycles in the data recorded by a spacecraft crossing the structure. Besides, (2) we anticipate that a scatter 
plot of V x versus the plasma density N during the passage of the whirling flow should show the statistical 
trend indicated in Figure 1 (bottom). Features (1) and (2) are the basic elements of a criterion that permits 
the identification of a boundary rollover in the observations. 

Aside from (i) a case study addressing an interval of southward pointing IMF [Hwang et al., 201 2a] and (ii) 
another study with a dawnward pointing IMF [ Hwang et al., 201 2b], most of the works on vortical struc- 
tures at the magnetopause/boundary layer have concentrated on a strongly northward pointing IMF, which 
is parallel to the Earth's field at low latitudes. If this lasts for several hours, it is typically associated with the 
northward pointing phase of interplanetary magnetic clouds [Burlaga etal., 1981]. A northward orientation 
favors the development of KFH instability because when the wave vector k of the perturbation is orthogo- 
nal to the average direction of the two magnetic fields ("flute" modes), or normal to the stronger one, the 
restraining magnetic forces are nearly canceled. At the same time, a substantial part of the velocity shear 
effect is retained. This argument applies equally well to southward IMF. 

Flere, by contrast, we focus on a situation where the IMF is oriented in a radial direction pointing approx- 
imately sunward, opposite to the solar wind flow. Under this configuration we present an example of 



-2 -1 1 


Figure 1 . A schematic to help interpret the Wind data, (top) The wavy 
magnetopause at the equatorial dawn flank is shown as it begins to roll 
over into a vortex by the KH instability. The magnetosheath flow is tail- 
ward (V x < 0), while the magnetosphere is stagnant. The drawing is 
shown in the frame of the vortex so that the cold dense magnetosheath 
tongue (blue) protruding to the left is slowing down relative to the aver- 
age flow, while the related hot tenuous magnetosphere is moving faster, 
(bottom) Expectations drawn from this for the scatterplot of V x versus N. 
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rolled-up flow vortices making up the entire boundary layer at low latitudes a few R E (Earth radii) tail- 
ward of the dawn terminator. The criterion for inferring the rolling-up stage, which was mentioned above, 
is satisfied. Furthermore, the Wind probe that recorded the rolling motion on 24 October 2001 was trav- 
eling orthogonal to the bulk motion of these structures, an ideal circumstance and one which is much 
superior to magnetopause-skimming orbits, which do not sample the whole structure of the vortices. In 
addition, the external field, too, was exceptionally steady and smooth in a plasma of low beta. In particular, 
there were no significant variations in the solar wind dynamic pressure. As noted by Farrugia etal. [2007] 
the magnetosphere was in a very quiescent state. Reconnection processes were at best weak and patchy 
(in time). 

We then inquire into the possibility that the vortices are of KH origin. We adopt two approaches. In the 
first approach, we input parameters to the theoretical stability analysis taken from the exact IVIFiD solution 
derived by Spreiter and Rizzi [1 974] and appropriate for collinear field and flow. This theory was applied to 
the present event during the later time when Wind was crossing the magnetosheath [Farrugia etal., 2010]. In 
the second approach we input to the theoretical calculations the observations made by the Wind spacecraft. 
In both cases we work with compressible MFID equations, using for the physical quantities continuous pro- 
files across a thick boundary layer. This avoids pitfalls in the use of the stability condition for a thin boundary 
model (Appendix A, formula (A1 2)), pointed out by Gratton et al. [2004a]. In both cases we find the region to 
be KFH unstable. This is thus the first reported instance of rolled-up KH vortices populating a boundary layer 
under a near-radial IMF. 

A magnetic field aligned with the flow is, of course, prima facie the most unfavorable configuration for 
the Kelvin-Helmholtz instability because the magnetic tension exerts a stabilizing action, which cannot be 
avoided by modes of the flute type. The stabilizing action of the field tension is precisely that avoided by 
flute modes. With a wave vector k perpendicular to the magnetic field, flute modes eliminate its operation, 
but in field-aligned flows they eliminate also the instability driver. In this paper we discuss how, nonetheless, 
the configuration can be KH unstable. 

A distinctive aspect of the case we present is that a radial magnetic field is forced to be drawn along by the 
billows when they arise. This constitutes a substantial difference from the KH instability for northward point- 
ing fields, where vortices can grow in a flute mode configuration, with only small changes in the orientation 
of the field lines. 

The layout of the paper is as follows. After discussing the interplanetary data, we describe the observations 
in the boundary layer made by Wind. We then discuss elements of the KH instability relevant to our work. A 
summary and discussion follows. We give some technical details on the KH instability in the two appendices. 

2. Observations 

2.1. Interplanetary: ACE 

Interplanetary conditions during the period we study consisted of a structure which formed the last in a 
set of interacting interplanetary coronal mass ejections (ICMEs). The ICMEs passed Earth during the 5 day 
period from 21 October to 25 October 2001 [see Farrugia et al., 2007] and were preceded by a strong shock. 
The state of the magnetosphere went from being strongly disturbed (21 October to 23 October) to being 
almost quiescent. In the first period, two intense geomagnetic storms (Dsf < -150 nT) were recorded. Then, 
on 24-25 October, all organized activity subsided. This very quiet period ended when a trailing shock was 
seen advancing into the ICME sequence. 

Observations over that part of this interval which is relevant to our study are shown in Figure 2. The 
interplanetary plasma and magnetic field observations are from the ACE spacecraft in orbit around 
the LI Lagrangian point. They were acquired by the Solar Wind Electron, Proton, and Alpha Monitor 
[McComas et al., 1998] and Magnetic Fields Experiment [Smith et al., 1998] instruments, and are at 64 
s (plasma) and 1 6 s (magnetic field) temporal resolution. The time interval shown is 1 800-2100 UT, 

24 October 2001 . From top to bottom the panels display the proton density, temperature, (in red: the 
expected temperature after the statistical analysis of Lopez [1987]), bulk speed, the GSM components 
of the magnetic field (color coded), the total field strength, the IMF cone angle, i.e., the angle made 
by the magnetic field to the Earth-Sun line, the dynamic pressure, the angle (shear) between the field 
and flow vectors, the proton beta, and the sonic and Alfven Mach numbers (M s and MA, respectively). 
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18 19 20 21 

UT 


Figure 2. ACE plasma and field observations during 1 800-21 00 UT, 24 October 2001 . The panels show the proton density, temperature (in red: the expected 
temperature for normal solar wind expansion), bulk speed, the total field and (colored) its GSM components, the IMF cone angle, the dynamic pressure based on 
the protons, the angle between the field and flow vectors, the proton beta, and the sonic and Alfven Mach numbers. 
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WIND 19UT, OCT 24-02 UT, OCT 25, 2001 (GSM) 



X(R e ) 


This period is marked by steady condi- 
tions and very smooth field and plasma 
temporal profiles. The temporal varia- 
tions, which were a leading feature of 
the previous 3 days, have died down 
completely. The data show a slow (aver- 
age and standard deviation: <V> = 
372.5 ± 2.5 km s _1 ) and very cold (<T> 

= 41 87 ± 375 K) ICME, the proton tem- 
perature being about 8 times less than 
the expected temperature. Compared 
to the normally dense slow solar wind, 
the density (<N> = 3.74 ± 0.46 cm -3 ) 
is about one half of a typical value of 
7-10 cm -3 , leading to below-average 
dynamic pressure of 1.0 ± 0.10 nPa. 

As a consequence of this, the proton 
/? is also very low, whence the smooth 
magnetic field profile. Because of the 
low T p , the sonic Mach number M s is 
very high, of order 35. At the position 
of Wind tailward of the dawn termina- 
tor, we therefore expect that effects due 
to compressibility of the plasma will be 
accentuated. The Alfven Mach num- 
ber is not particularly small ( — 7), and so 
the magnetic forces should not have a 
dominating influence on the instabil- 
ity. Importantly, the magnetic field has a 
near-radial orientation (Figure 4, fourth 
row) with a cone angle of 17.7° ± 4.1°. 

It makes an angle with the plasma flow 
vector of 1 62.0° ± 4.7° so that it points 
almost opposite to the solar wind. 

2.2. Observations in the Dawnside 
Boundary Layer: Wind 

In October 2001 the Wind spacecraft 
was orbiting the magnetosphere, reach- 
ing perigee in the near-geomagnetic 
tail region. Figure 3 shows its orbit 
from 1 900 UT (24 October) to 0200 UT 
(25 October), after which time it exited 
into the solar wind [Farrugia eta!., 2010]. 
The red segments indicate the time 
when Wind was traversing the dawnside 

low-latitude boundary layer (LLBL) (1900-20:30 UT, see below) downstream of the terminator at X 13.5 

R e and at somewhat northerly GSM latitudes (Z ~ 5.5 R E ). This orbit cuts across any structures which are 
propagating downstream in this region. This is an ideal situation for our purposes. 


N 



X(R e ) 

Figure 3. A segment of Wind's orbit for the time interval 1 900 UT 
(24 October)-0200 UT (25 October). The plot shows the trajectory 
projected into the X-Y and X-Z planes. Tick marks are shown at each 
hour. The red segment (19:00-21:30 UT) refers to the time Wind 
was crossing the low-latitude boundary layer (LLBL), thus moving 
predominantly dawnward. 


Examples of the structures encountered in the period 19:00-19:30 UT are shown in Figure 4. The data are 
from the 3-D Plasma analyzer (3DP) [Lin et al., 1 995] and the magnetic field investigation [Lepping et a!., 
1995], both plotted at 3 s resolution. Shown from Figure 4 (first-tenth rows) is the proton density, bulk 
speed, temperature, the total field and its GSM components, and the GSM components of the flow vector. 
The dashed blue line in Figure 4 (second row) gives the average magnetosheath velocity in the first half 
hour after Wind's entry at 20:30 UT (= 314 km s -1 ). The averages of the bulk flow velocity components are 
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WIND/MFI/3DPL October 24, 2001 19:00-19:30 UT (GSM) 



19: 00 05 10 15 20 25 30 


UT 

Figure 4. Proton plasma and magnetic field observations from Wind for the period 1 9:00-1 9:30 UT. (first-tenth rows) the proton density, bulk speed, temperature, 
the total field and its GSM components, and the GSM velocity components. The dashed blue line in Figure 4 (second row) gives the average magnetosheath 
speed. Dashed red lines in Figure 4 (eighth-tenth rows) show the average values of the respective quantities over the interval plotted. Note the speeds of the hot 
tenuous plasma, which exceed the solar wind. 
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marked in Figure 4 (eighth-tenth rows) by the horizontal red lines. The average negative V y component 
(= -46 km s _1 ) results from the dawnward flaring of the magnetopause (« 7.5°) at Wind's dawnside locale. 
Marked by vertical guidelines are the times when sharp increases in density occur and when simultane- 
ously impulsive changes in the magnetic field and/or plasma parameters are evident. It should be noted 
that uncertainties in the plasma moment results do increase as the proton density drops below ~0.1 cirr 3 . 
However, nowhere are the interpretation and conclusions of the paper affected by this. 

We note the following features: 

1. Repetitive high-speed bursts of a hot tenuous plasma reaching speeds (up to ~650 km s -1 ) which are well 
in excess of the magnetosheath speed. 

2. After the discontinuities (vertical lines), intervals of a cold dense (magnetosheath) plasma each lasting for 
~2-3 min are encountered. 

3. In the Earth's frame, the cold dense plasma is moving more slowly tailward than the average flow. In the 
average velocity frame, its motion is thus predominantly sunward. 

4. By contrast, in the average velocity frame the hot tenuous plasma is moving antisunward. Note the 
repeated overshoot of this plasma with respect to the antisunward velocity. 

These last two points may be seen very well from the clear anticorrelated behavior of the density N and the 
antisunward velocity, —V x . 

5. Sharp changes in the field and flow vectors, including abrupt polarity reversals, tend to occur at the 
leading edges of the cold dense (magnetosheath) structures. With one exception (that at 1 9:1 1 UT) the 
leading edges are thus simultaneously current and vortex sheets. 

6. Considering only the largest changes (indicated by the vertical guidelines), there are six intervals of 
roughly repetitive structures with an average duration of ~4.5 min. This average periodicity is retained 
throughout the entire 1 .5 h traversal of the LLBL, as we discuss below (section 2.3). 

From (1) to (4) we see that the criterion for identifying rolled-up vortices given in the section 1, namely, a hot 
tenuous plasma flowing at speeds higher than that of the magnetosheath, is well satisfied. We exclude the 
possibility that the faster-than-sheath flows of the hot tenuous plasma resulted from reconnection at the 
tail flank magnetopause. The Wind location in the Northern Hemisphere tailward of the terminator gives a 
generally sunward pointing field (B x >0) which is approximately parallel to the interplanetary field. So the 
large magnetic shear required for reconnection to occur at the tail flank magnetopause near the dawn-dusk 
terminator [Gosling eta!., 1986] is not present. 

We now illustrate these features by focusing on one typical cycle of the plasma and field behavior. Plasma 
and magnetic field data for a single cycle, corresponding to the interval marked by the horizontal red 
bar in Figure 4 (first row), are shown in Figure 5. The same quantities as in Figure 4 are plotted in Figure 5 
(first-seventh rows; note, however, the linear scale for the density). Figure 5 (eighth-tenth rows) shows the 
plasma velocity in the average velocity frame. The approximate duration of the cycle is from — 1 9:1 2:20 UT to 
— 1 9:1 8:00 UT (~ 5.7 min). The cold dense plasma interval is bracketed by the two vertical dashed red lines 
and lasts for ~2.8 min. Immediately preceding the leading edge of the cold dense plasma at 19:12:20 UT, a 
plasma of low density and elevated temperature is moving at a speed exceeding that of the solar wind. The 
plasma there is flowing mainly perpendicular to the local magnetic field (not shown). The same may be seen 
from ~1 9:1 6 to 19:18 UT ahead of the next cold dense plasma burst. 

Relative to the average velocity, the cold dense plasma is moving mainly sunward (Al/ X > 0). With an aver- 
age speed of 280 km s -1 and a duration of 2.8 min, the scale size of the cold dense plasma in the X direction 
is estimated as 7.4 R E . So it is very stretched in the X direction, compared, say, to the distance around the 
magnetopause from the nose which is of order 20-25 R E . After the cold dense plasma there follows a stage 
(~ 1 9:1 5:1 0—1 9:1 6:00 UT) in which the plasma has acquired a dawnward velocity component (AV y < 0) 
and its sunward speed has decreased, then comes a burst of hot tenuous plasma moving antisunward and 
northward and which ends up moving strongly antisunward and duskward. A rotational motion superposed 
on the antisunward and dawnward bulk flow is thus evident. 

To visualize the flow rotation in the average velocity frame, we show in Figure 6 the residual flow vectors 
AV X , AV y for the period 1 9:1 2:30-1 9:1 8 UT. Time runs from the bottom to the top, and the labels "S" and "E" 
refer to the start and end of the interval. The arrows show the coordinates of the residual vectors. The blue 
arrows refer to the cold dense plasma, the red arrows to the hot tenuous plasma, and the green arrows to 


FARRUGIA ETAL. 


©2014. American Geophysical Union. All Rights Reserved. 


4578 


IgAGU Journal of Geophysical Research: Space Physics 


10.1 002/201 3JA01 9578 


WIND/MFI/3 DP October 24, 2001 19:12:00-19:18 UT (GSM) 



19 : 12 13 14 15 16 17 18 

UT 

Figure 5. Wind plasma and field data for the time interval ~1 9:1 2 UT to ~1 9:1 8 UT. The format is the same as for Figure 4 except that the last three rows show the 
velocity components in the average velocity frame, i.e., when the average velocity computed over this interval is subtracted. 
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an intermediate state in ( N , T). It is seen 
that, in the average velocity frame, the 
cold dense plasma is flowing mainly 
sunward with generally a very small 
duskward component (AV y >0). The 
hot tenuous plasma (red) first flows tail- 
ward and then tailward and duskward. 
The flow direction of the plasma in 
between (green) starts rotating from 
a sunward and dawnward orientation 
and finishes in an antisunward orien- 
tation. This provides clear evidence of 
rolling up (see section 1). 

2.3. General Features of the 
Vortical Structures 

The fact that the hot tenuous magne- 
tospheric plasma is moving at speeds 
above those of the solar wind is strong 
evidence that the structure we are 
dealing with in Figure 6 is a rolled-up 
vortex [Nakamura et a!., 2004; Takagi et 
al., 2006; Hasegawa et a!., 2004, 2006]. 
We wish now to confirm this for all 
the quasi-periodic structures seen by 
Wind in the interval 19:00-20:30 UT. 

In Figure 7 (bottom) we plot V x versus N for the whole interval. Figure 7 (top) shows the same quantities 
for the one plasma and field cycle we have just been discussing. The color scheme indicates the tempera- 
ture (red for hot and blue for cold). The dashed horizontal lines show the corresponding quantities when 
Wind entered the magnetosheath, (not shown) [see Farrugia et al., 2010, Figure 7], averaged over the first 
half hour. 

From Figure 7 (bottom) it is seen that the bulk of the hot tenuous plasma is moving in an antisunward direc- 
tion faster than the magnetosheath. Clearly, also, the figure shows that the origin of the cold dense plasma 
is the magnetosheath. This is the same trend as seen in the single cycle plotted in Figure 7 (top) although 
the highest speeds recorded there were ~ 580 km s -1 . Following Takagi et al. [2006], we conclude that Wind 
is observing an LLBL populated entirely by a sequence of rolled-up vortices. 

Figure 8 depicts the motion of the plasma in the dawn-dusk direction in the form of a scatterplot of the 
residual A1/ y versus AV X for the whole interval 1 9:00-20:30 UT. The color is proportional to log T (red = hot) 
and the size of the squares to N. The figure shows a continuous distribution of Al/ y values spanning across 
zero. There is no strong preference for positive or negative Al/ y .The spread in AV y of the hot tenuous plasma 
is wider. This overall picture confirms the persistence of the rotational motion in the average velocity frame 
quite clearly. 

We recall from Figure 2 that interplanetary conditions were steady. Specifically, there were no significant 
variations in the dynamic pressure, P dyn . But, in fact, there are large-amplitude, quasi-periodic fluctuations 
of this quantity generated by the vortices themselves. To show this, we consider in Figure 9 the temporal 
variation of P dyn at Wind. The 1.5 h traversal is split into three ~0.5 h segments which are plotted underneath 
each another. Each panel shows the thermal plasma pressure (green trace) and the dynamic pressure (black 
trace). The blue and red traces are 21 point (~ 1 min) running averages of these two quantities, respectively. 
In computing the dynamic pressure, we took into account the average flaring of the dawn magnetopause, 
given by V y in Figure 4. This gives a flaring angle of 0 « 7.5°, and the dynamic pressure has been multiplied 
by sin(P). 

One can see that the average dynamic pressure is subject to large-amplitude oscillations of period ~5.0 min. 
Note that there are six clear waves corresponding to the vortices in the Figure 9 (top). These are the 
ones identified in Figure 4 except for the small one between 1 9:1 1 and 1 9:1 2.2 UT, which might indicate 



Figure 6. Residual vectors in the XY plane for the structure in Figure 5. 
Symbols S and E mark the start and end of the structure. The labels CD 
and FIT refer to "cold dense" and "hot tenuous," respectively. Time runs 
from bottom to top. The blue, green, and red vectors represent differing 
plasma parameters, as explained in the text. In the average velocity frame 
shown, the flows start moving sunward and slightly duskward (blue). 
They then rotate dawnward and become progressively antisunward 
(green) and finish flowing antisunward and duskward (red). 
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Figure 7. Scatterplots of (top) V x versus N for the vortex at 19:12:30-19:18:00 UT and (bottom) for the whole LLBL 
crossing. The logarithm of the temperature is indicated by the colors, where red = hot and blue = cold. Velocities are 
plotted in the Earth's frame. The horizontal dashed line marks the average magnetosheath speed observed when Wind 
crossed into this region at 20:30 UT (not shown) [see Farrugia et at., 201 0, Figure 7]. Figure 7 (top and bottom) shows 
the presence of (i) a hot tenuous plasma moving at high speeds tailward and (ii) a dense cold plasma moving antisun- 
ward at speeds close to that of the magnetosheath. The figure shows clearly that the origin of the cold dense plasma is 
the magnetosheath. 


some ongoing coalescence. The thermal plasma pressure behaves as the dynamic pressure, only at much 
reduced amplitude. 


Underneath each plasma pressure panel, we plot the magnetic pressure (P 8 ) for the corresponding interval. 
A linear scale is used for this. Overlaid on the P 8 trace we show the total pressure (thermal + magnetic; P tot 
in blue). The curves practically coincide because the magnetic pressure is much larger than the thermal 
pressure. Overall pressure balance is not maintained. As Mura [1 997] pointed out, the quasi-equilibrium in a 
nonlinear vortex can be achieved if the centrifugal force exerted by the vortex flow is balanced by the total 
pressure gradient force. From this, an estimate of the total pressure variation in a vortex was obtained by 
Hasegawa [201 2]. It is 


AP tot > 0.5p(Al/) 2 
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LLBL Wind data, October 24, 2001 , 19:00 - 20:30 UT 


where p is the mass density, and AV 
is the velocity jump. Estimates of 
these quantities in our case are N = 
0.3-0.5 cm -3 and AV = 300 km s _1 . 
We thus require that AP tot to be of 
order, or larger than, 0.02-0.04 nPa. 
This is roughly also what we have 
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(see Figure 9). 

Lastly, we note that the fluctuations of 
the thermal pressure can produce ion 
acoustic waves along the geomagnetic 
field. The variation of the magnetic 
pressure can radiate magnetosonic 
waves across the magnetic field. The 
vortices can thus give rise to large-scale 


-400 -300 -200 -100 0 

A V , km/s 

X 


100 200 


300 effects in the plasma sheet. 


Figure 8. The figure shows a scatterplot of Al/ y against AV X for the whole 
interval. Color is proportional to log T, and size is proportional to N. The 
plot is in the average velocity frame. The distribution in the dawn-dusk 
direction AV y shows a wide spread across zero. It is wider for the hot 
tenuous plasma. 


3. Generating Mechanism: 
The Kelvin-Helmholtz 
Instability Source 

The configuration of 24 October 2001 
appears not to favor the onset of the 


KH instability: (i) the restraining magnetic forces are strong in field-aligned flows and (ii) the large M s ushers 
in the other stabilizing factor, compressibility. 

We now examine the issue more closely. We model the LLBL transition by continuous functions for the phys- 
ical parameters. We call this "thick model" for short, to distinguish it from a "thin" approximation where the 
quantities suffer a discontinuous change across the boundary layer. (Appendix A, (A1 2)) [see also Grafton et 
al., 2004a, 2004b; Gnavi et a!., 2009]. For the stability analysis we work in a flow-aligned coordinate system 
defined as follows. The x-axis points in the direction of the local V. They-axis points across the LLBL, normal 
to the local magnetopause and directed outward. The z axis completes the right-handed Cartesian triad and 
is oriented in the same sense as geomagnetic north. 

Scalar and vector quantities in the LLBL are represented by hyperbolic tangent functions with a scale length 
d, for example, 


for the velocity, and with similar expressions for 8 and N (Appendix A). Subscripts "1 " and "2" refer to mag- 
netosheath and magnetosphere quantities, respectively. The temperature profile T(y) follows from the 
pressure balance equation across the layer (Appendix B). We take D = 4d as a representative value for the 
LLBL thickness, which ranges approximately from y = -2d to +2d. The normalized quantities contain d and 
1/, implicitly such as, for example, in the normalized growth rate g = yd/1/,. An estimated value of D, and a 
measured value of V v can be introduced in the discussion at the end of the theoretical calculation; it is not 
necessary to assume them beforehand. The compressible MHD stability theory used here is summarized in 
Appendix A, where some details of the procedure can be found. For every k mode the KFH instability is driven 
by the intensity of the velocity projection V k in the k direction (V k = V -k, where k = k /\ k |). The mag- 
netic tension that opposes the instability depends on the magnetic field projection B k in the k direction 


We follow two approaches regarding the physical parameters which we input into the theory. In the first 
approach, the Mach numbers are based on the Spreiter and Rizzi [1 974] theory that gives an approximate 
representation of the solar wind-magnetosheath transition for collinear MFID flows. We use the Spreiter - 
Rizzi solution with solar wind input from ACE. Close to the terminator this theory predicts approximately the 
following: M s = 7.7, MA = 4.9, which corresponds to a magnetosheath plasma /?, = 0.97. 


l/ x = i/,(i +tanh(y/d))/2, 


(D 


(B k = 8 -k). 
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Figure 9. For three subintervals, the plot shows pairwise the proton thermal pressure and dynamic pressure and below the magnetic pressure and the total 
pressure (blue trace). The smoothed averages of the thermal and dynamic pressure are shown by blue and red traces, respectively. 
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The geomagnetic field on the dayside 
is assumed to be perpendicular to V , 
i.e., the magnetic shear angle in this 
model is 90°. Presumably, this angle 
was not exactly 90°, and it varies with 
the distance from the subsolar point. 
But we think that near-Earth deviations 
from 90° could not have been substan- 
tial. Anyway, the chosen shear angle is 
not critical to decide on the instability 
because we intend to switch off — or, 
at least, much reduce — the magnetic 
tensions on the magnetosphere side of 
the LLBL by considering k vectors nor- 
mal to the local geomagnetic field. This 
choice of k favors V k , the driver of the 
instability (and maximizes it when the 
magnetic shear angle is exactly 90°) 
but exposes the k mode to the full sta- 
bilizing influence of the magnetosheath field projection B k . For the dayside we assume a typical particle 
density ratio A/j/A/, ~ 0.1. 

The pressure balance equation imposes an upper limit on the magnetic field ratio Bj/B, < 1 .4. (About this 
requirement, see condition (B7) in Appendix B.) In approach (1) we computed with B 2 /B , = 1, n 2 / n i = 0.1, 
and (as a consequence of equation (B6), Appendix B) a temperature ratio T 1 /T ] ~ 10. The choice reflects 
expected values at the MP away from the subsolar point, but still near Earth, as the terminator. 

The mode considered is with k parallel to the flow. (Computation shows it to be the k orientation of fastest 
growth.) Figure 10 shows the (normalized) imaginary part of the characteristic value, c, = y/f/d/,), as a func- 
tion of kd. Quantity c,- is a linear function of kd in most of the interval (for kd > 0.15) so that the growth rate 
g = yd/V-f as a function of kd is approximated by a parabola (not shown). The maximum of the normalized 
growth rate g = c,kd = yd/ 1/, occurs at kd =0.245 (7 = 6.41D) for g = 0.0832. From this value we may esti- 
mate an e-folding time r e = \/y ~ 48 s for LLBL sites near the terminator, with D = 0.5 R E , assuming that 
the LLBL thickness is not yet broadened by the instability, and with 1/, ~ 200 km/s. Therefore, in a bound- 
ary layer not yet widened by perturbations, the KFH instability can grow quite fast. We think that the vortices 
observed by Wind are generated in the LLBL closer to Earth. 

The values chosen for D and 1/, are not experimentally verified, because of the absence of in situ spacecraft 
measurements. Flere we could adjust to a certain extent the estimates, since they depend on assumptions 
of D and 1/, . In short, the values quoted are only indicative. 

We may note that the length of the new generated perturbations increases during the transport toward the 
tail. This is because of the difference of magnetosheath speed between the front part and the tail part of the 
(rather long) structures. Thus, the structures are stretched during their journey to the Wind locale. 

In addition, more changes of the vortex aspect ratio take place during the rollover, which modify the local 

boundary layer shape, and thickness. Important nonlinear strains of 
the structure occur, particularly when the radial solar wind field is 
affected and is dragged in, by the vortex. At this stage of the evo- 
lution, the physics indicates that structures with long scale lengths 
are favored, because of the smaller curvature of the magnetic lines, 
and the consequent lower magnetic tension that operates against 
the process. 

The two mentioned effects may account for a significant incre- 
ment of the wavelength observed at the Wind locale. Other possible 
causes of the long scale length structures observed in the flanks are 
discussed in Hasegawa [2012]. 


Table 1 . Input Parameters for the 
Stability Analysis 3 

Magnetosheath Magnetosphere 

M s = 5.6, MA = 6.8 
1/, = 314 km/s V 2 =0 

A/, = 5 cm- 3 N 2 = 0.05 cm" 3 

B, = 4.5 nT B 2 = 9.5 nT 

Tl = 180° Z 2 = 1 09° 

a The angle that the magnetic field 
makes with the x axis is denoted by /. 


October 24, 2001 , Thick Model 



Figure 10. Normalized imaginary part of c, i.e., the complex phase 
velocity y/(kVi), versus kd. Flyperbolic tangent model with input param- 
eters from the Spreiter-Rizzi theory (1974) for a boundary layer at the 
terminator. Quantities M s = 7.7, MA = 4.9. 
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October 24, 2001, Thick Model 



Figure 11 . Normalized growth rate as a function of kd. The input param- 
eters for the stability calculations are based on Wind data. Maximum 
growth rate is reached at kd = 0.43. For further details, see text. 


We now discuss the second approach. 
Here we input to the model data 
acquired during the Wind's LLBL traver- 
sal. The scenario is a composition of 
averages of measurements: (a) made in 
the magnetosphere before, but close 
to, the entrance to the LLBL (including 
early times of the passage through it) 
with (b) data recorded — albeit later — in 
the adjacent magnetosheath. 

At Wind's position, with a magne- 
tosheath average velocity of 314 km s -1 , 
the Mach numbers computed with 
data for that region are M s = 5.6 and 
MA = 6.8. They lead to a plasma of 
/?, = 3.5. The magnetic field B in the 
magnetosheath is still approximately 
collinear with the flow [Farrugia et at., 
201 0]. A difference from approach 1 is 
that the magnetic shear angle is not 
90° (as hypothesized for near-Earth 
positions) but 71° with respect to 
the sunward direction (taken from an 
average B on the magnetosphere side). 


To sum up, the input parameters for the stability analysis are presented in Table 1 . Inside the magnetosphere 
the total pressure is almost purely magnetic, but in the magnetosheath thermal and magnetic pressures are 
of the same order. The Bj/B, ratio satisfies the condition (20) given in Appendix B. The ratio Bj/B, comes 
close to the upper limit of that formula but stays within the validity range of the pressure balance condition 
across the boundary layer (BL). (Note also that the upper bound of condition (20) could be a bit larger than 
the quoted value, if a correction for flaring is accounted for). 


Figure 1 1 shows the normalized growth rate g a function of kd. It reaches a still significant maximum value 
g = 0.033 at kd = 0.43 and goes to zero at kd « 0.07 (long A) and kd = 0.84 (short A). The angle of k with 
the x axis is 4> =19°. With a ratio N 2 /N , = 10~ 2 the growth rate y is zero in the kd -> 0 limit. This is a case of 
stability at long wavelengths. 


The second approach intends to show that a steady state LLBL model, endowed with equilibrium quanti- 
ties represented by continuous functions that connect mean values of magnetospheric and magnetosheath 
data, is unstable. The averages include long stretches of time on either side of the LLBL, because the instabil- 
ity is found with wavelengths of several R E , and the penetration depth of the KH perturbation is expected to 
be large. However, the LLBL and the magnetosheath are both perturbed already. The former by the passage 
of vortical structures as discussed in section 2, which we conclude are formed at some place near Earth; the 
latter by the turbulence after the bow shock, and by large amplitude oscillations of long period [Farrugia et 
al., 2010]. Nonetheless the result is an indication that the LLBL, at the Wind's orbit position, has amplifying 
properties regarding the KH instability mechanism. In other words, if by reason of intermittency the pas- 
sage of vortical perturbations is temporarily suspended, the (unstable) background LLBL still maintains the 
capability to grow perturbations, and eventually to roll over the velocity gradient layer again. 


4. Summary and Discussion 

We took advantage of a rare coincidence of a long interval of radial IMF, steady solar wind conditions, and 
a spacecraft taking observations along a path that cuts perpendicularly through the near-Earth flank of 
the magnetosphere. The Wind observations through the LLBL at X = -13 R E showed the LLBL to be full of 
rolled-up vortices. These were shown to arise from the KH instability. The new result here is not so much the 
observation of KH-like oscillations but that they occur under a radial IMF, which should suppress the growth 
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of KH waves. So while rolled-up vortices for northward IMF have been reported before, this is first reported 
case of observations of KH rolled-up vortices for radial IMF. 

In some ways this was a continuation of work started in Farrugia et al. [2007, 2010], In those papers we 
focussed on the extremely quiet state of the magnetosphere after a 3 day long period of strong disturbance. 
The cause was a series of ICMEs, and the period we concentrated on here constituted the last of these where 
very steady conditions prevailed. The Farrugia etal. [2010] study of this event dealt with the entire magne- 
tosheath showing that the near-parallel alignment of field and flow held throughout the magnetosheath 
and that the data matched a relevant theory that treats flow-aligned field in the magnetosheath. 

In a 1 .5 h traversal we identified approximately 1 5 vortices. We argued they had reached the nonlinear 
stage and had started to roll up. The rolling-up process was inferred from the repeated presence of a 
low-density, magnetosphere plasma moving antisunward at speeds greater than in the magnetosheath, 
which recent studies have shown to be a reliable indicator of such structures based on single-spacecraft in 
situ observations. 

We then presented two KH instability calculations using different inputs for the theory model of the tran- 
sitions. In the first we inputted data from the theory of Spreiter and Rizzi [1974], which gives an exact MHD 
solution for field-aligned flows. In the second, we inputted direct measurements made by Wind in the mag- 
netosphere and in the magnetosphere at the beginning and end of the LLBL crossing, respectively. In both 
cases the LLBL was found to be unstable. 

Although the solar wind dynamic pressure was very steady, the passage of the structures gave rise to 
large-amplitude modulations of the magnetic pressure and the dynamic pressure in the boundary layer. 
This could set up waves traveling along (ion acoustic waves) and perpendicular (magnetosonic waves) to 
the magnetic field. This shows that the passage of the large vortices at dawn could influence large parts of 
the plasma sheet. 

We now discuss various aspects of the observations of the 24 October event. In previous work, Farrugia et 
al. [2007] concentrated on the very low level of geomagnetic disturbances which prevailed on this day. Such 
were, for example, an average Kp index = 0+, polar caps which had very weak electron precipitation without 
any consistent north-south asymmetries, and patchy and weak reconnection at low latitudes or poleward 
of the cusp. In particular, the authors noted a cross-polar cap potential (CPCP) of ~20 kV. This would be 
of about the same magnitude as that commonly ascribed to the contribution to the CPCP of viscous-type 
interactions. However, Farrugia et al. [2007] argued against the KH instability. This conclusion was essen- 
tially based on the lack of ULF pulsations of the geomagnetic field in the Pc5 range (2-7 mHz), which the 
KH instability is often thought to give rise to via the field line resonance theory [Chen and Fiasegawa, 1974; 
Southwood, 1 974], They thus could not find a solar wind driver for the weak and patchy convection. 

Was the instability a result of the (radial) direction of the IMF, or was it favored by specific values of the Mach 
numbers in this case? This is an important question that deserves further attention. We think that in this case 
MA values (estimated at the terminator and measured at Wind's orbit) were clearly helping the onset, and 
development of the instability. The radial field orientation is unfavorable, in general, to the KH instability. 
However, even in the case of normal solar wind and Parker's spiral field, there is the possibility that KH waves 
may develop inside the boundary layer, and grow thereafter tailward. Further work is necessary to test other 
cases to see if it was the radial IMF orientation that did this, or rather other favorable parameters. 

From the theory it was concluded that (i) in both approaches the boundary of this wide LLBL was KH unsta- 
ble and (ii) the long wavelength limit is stable. That is, a thin layer would be stable. The instability appears 
only with the thick boundary layer. We also concluded that the generating site was well upstream of the 
observation locale. As a consequence of the stability study, we assume, as seems reasonable, that the life- 
time of each member of the vortex sequence is similar because they are generated approximately at the 
same position upstream (closer to the Earth). 

If a magnetic field collinear with the flow is generally unfavorable to the development of the KH instability, 
and the magnetosheath flow was supersonic, why is the LLBL unstable in this event? This is so because the 
physical conditions of the inner edge of the LLBL are very different from those of the outer edge (adjacent 
to the magnetosheath). That the inner edge of the LLBL, in general, may be prone to the KH instability was 
first pointed out by Sonnerup [1 980] and has been considered by other authors, among them Ogilvie and 
Fitzenreiter [1 989], and Miura [1 992]. 
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Behind the bow shock, the magnetosheath is frequently in a turbulent state. The plasma of the boundary 
layer is pulled along by the solar wind. A velocity shear flow parallel to that of the magnetosheath is estab- 
lished across this layer. The motion is subsonic inside the LLBL, because of the decreased speed with respect 
to the magnetosheath, and the higher temperature of the magnetospheric plasma. At the equatorial day- 
side, the geomagnetic field is mainly normal to the flow. At the inner edge side, even if moving at reduced 
speed, the obstacles to the growth of the KH instability are attenuated. Flute modes with a wave vector nor- 
mal to the local geomagnetic field, and parallel to the internal flow direction, are not restrained by magnetic 
tension forces, and the low compressibility reduces additional stabilizing effects. Downstream, MA increases, 
and magnetic tensions are further reduced so that conditions for instability improve. 

A major reason why people are interested in the nonlinear stage of the KFH instability in the first place is that, 
by breaking the frozen-in condition, it offers the possibility for mass transfer. This transfer would happen at 
current sheets where oppositely directed magnetic fields have been brought next to each other during the 
rolling-up process (see Figure 4). We found several instances of current and vortex sheets. Indeed, most field 
and flow cycles contained one of these. 

Current sheets are prone to magnetic line tearing and hence are a possible way to mass transport across the 
MP (see, e.g., Otto and Fairfield [2000] and Otto and Nykyri [2003] for computer simulation studies of field 
lines coiled up inside vortices). From the data we cannot tell to what extent the field lines are entrained by 
the vortices. But we find repeated evidence of current sheet formation. Whether mass transport is actually 
taking place in our case will be pursued in a further study. 

Our work and that of Hwang et al. [201 2a, 2012b] show that it is not necessary to have a northward point- 
ing IMF to excite the KFH instability. Neither is it necessary for the IMF to point north to produce field 
configurations conducive to reconnection and, by implication, mass entry. 

Appendix A: KH Theory 

The LLBL model with continuous functions used in section 3 describes a MHD-parallel flow with a local x axis 
directed along the velocity field. (The flow does not change direction.) The physical quantities are constant 
over (x,z) planes and vary only in the transverse y direction, chosen to be normal to the MP. Flowever, in 
general, the magnetic field may change both in direction and strength. The unperturbed (or average state) 
LLBL model is given by a set of functions: V = (V x (y), 0, 0), B = (B x (y),0,B z (y)), p(y) = m p N(y), for velocity, 
magnetic fields, and mass density p or particle density N, respectively. The temperature function T(y) results 
from this set of functions and the pressure balance equation (Appendix B). 

Across the LLBL, the physical quantities have hyperbolic function profiles: 


V x = V,(1 +tanh(y/d))/2, (Al) 

B = (B, + B 2 )/2 + (S, - B 2 )tanh(y/d)/2 (A2) 

9 = (0, + 0 2 )/2 + (6>, - <9 2 )tanh(y/d)/2 (A3) 

B x — B cos(B), B Z = B sin(B) (A4) 

N = (A/t + N 2 )/2 + (A/, - /V 2 )tanh(y/d)/2 (A5) 


where d is a scale length. The width D of the thick LLBL model is taken as D = Ad. 

The perturbation modes of the KFH instability are of the form 

E = f(y) exp(-iwt + ik x x + ik z z), (A6) 

where H is the y component of the Lagrangian displacement of a plasma element from a steady state posi- 
tion and f (y) is the corresponding amplitude. The (complex) angular frequency of the modes is denoted 
by co = co r + iy. The real part co r gives the frequency of the oscillations, and the imaginary part is the 
growth rate of the instability when y > 0; the e-folding time is i e = 1 /y. The wave vector is represented by 
k = (k x , 0, k z )\ k = | k | is the wave number; A = 2%/k is the wavelength. 
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The amplitude of the Fourier modes of the KH perturbation is governed by the second-order 
differential equation, 


H( 1 - — ) ^ 
dy V M ) dy 


- k 2 Ht; = 0, 


(A7) 


derived from the linearized equations of ideal (nonresistive), compressible MHD [ Gratton et at., 1 988]. A com- 
plex phase velocity c = co/k is introduced so that the functions H(y) and M{y) of the differential equation for 
f can be written as 


H(y) = p[(c-V k ) 2 -V 2 Ak ], 


(A8) 


M = 1 - 


C s+ V A 

(c-v k y 


+ 


cl/; 


Ak 


(c-v k r 


(A9) 


where V A = 8/ yjA-np is the Alfven speed, V Ak = B k /\jAnp, is a projected Alfven speed, c s is the speed of 
sound, and V k , B k are projections of the velocity and magnetic fields in the k direction. All these quantities 
are functions of y. The analysis is of a temporal type, that is, a (real) wave number k is given (as a Fourier 
component of the initial perturbation), and the response of the system determines the unknown (complex) 
value of c. To obtain c a boundary value problem for equation (A7) must be solved. 


When c s -»• oo, the coefficient M -y oo, and equation (A7) reduces to 


d_ L d£ 

dy dy 


- k 2 H£ = 0. 


(A10) 


that represents the incompressible MFID approximation. When the transition layer is very thin with respect 
to the wavelength, that is, when kd <s: 1, an approximate dispersion relation can be derived: 


H, + H 2 = 0, 


(All) 


where H, and H 2 are the values taken by H on each side of the BL; labels 1 and 2 refer to the magnetosheath 
and magnetosphere, respectively. This is the"thin model" result for incompressible plasma flows. From 
equation (A1 1) a well-known stability condition follows, 


p R (Ai;.k) 2 <i-[(B;.k) 2 + (8;.k) 2 


(A12) 


where p R is defined by 1 / p B = (i/p, + 1 / p 2 ), and A V = I/, - V 2 . The thin model condition, often used in 
the current literature, ensures stability when it holds for all directions of k . The thin model stability does not 
depend on the wavelength. 

The intricacy of the boundary value problem for equation (A7) with finite wavelengths derives from the fact 
that c is not an eigenvalue but a characteristic value entangled in a nonlinear fashion in the functions H and 
M. Moreover, when the direction of k changes, the functions V k (y) and B k (y) (and other functions, such 
as c s ( y)) also change. Thus, the analysis requires the solution of separate differential equations for every k 
direction. In this paper we solved the boundary value problem for c using a conventional shooting method. 
The compact form of equation (A7) facilitates the use of shooting methods. 


Appendix B: Pressure Balance Condition 

The field functions of the local LLBL model, B(y),p(y), T(y), etc., must satisfy pressure balance, 

B 2 B(y) 2 B 2 

P, + ^=P(y ) + ^ r =P 2 + ^. (Bi) 

Flere p = nk B (Ti + T e ) is the thermal pressure (/c 8 is Boltzmann's constant). We assume a common temperature 
value T-, = T e = T (in our case the proton temperature is from spacecraft data). 

It is convenient to write equation (B1) in terms of M s and MA, both computed with magnetosheath parame- 
ters adjacent to the local LLBL. Quantity M s = l/,/c sl , where c sl = i \J(yk B T e /m j ) = JiykgT^ /m p ) {y = 5/3, 
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and m, = m p is the proton mass). Similarly, MA = with V M 

equation (B1) can be written in the form 


P i 


p(y) + 


B(y) 2 

8rr 


P^ 


-c](y)+V 2 .{y) 


B^/ -y/ 4icp 1 = 8, /-y/Ajtn, m p . Then 


fl-L + J-V 

V Y M] MA 2 ) 


(B2) 


When M s and MA are known, the plasma beta is fixed because 

2nk B T = 4f[ 

P B 2 /8n y V 2 A ’ 

and since we are interested in the magnetosheath beta, 

_ 4 MA 2 
y M]' 

From these, the temperature function across the boundary layer can be written as, 


(B3) 


(B4) 


ny) _ ”i L 1 / B(y) 2 V 

n(y) [ /?, y B 2 J 

The local magnetosphere-to-magnetosheath temperature ratio is therefore 


T 2 _ n. 

1 + 1 

MU 

T^ n 2 


»- }\ 


(B5) 


(B6) 


which implies that the magnetosheath beta, /?,, together with the magnetic field intensity ratio B 2 /B 2 , sets a 
limit to the steady state boundary layer models. A local pressure balance does exists when 


8 2 /B 2 <1+^, (B7) 

and we see that T 2 becomes zero when B 2 /B 2 = 1 + /?,. Under ordinary conditions, the magnetosheath/?, is 
much larger than unity so that this limitation is not important. But in the 24 October 2001 event the values 
of/?, are comparable to, or even smaller than, unity. Flence, when setting stability models the constraint (B7) 
must be taken into account. 

Condition (B7) needs a correction when the boundary is flared with respect to the solar wind flow due to the 
presence of a normal component of the momentum flux. In practice, this can be approximately assumed as 
an increment of the effective B 2 , and then (B7) becomes a less severe bound. 
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